The gradual destruction of magnetism in the superconducting family NaFei^Co^As 
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The interplay and coexistence of superconducting, magnetic and structural order parameters in 
NaFei-^Co^ As has been studied using SQUID magnetometry, muon-spin rotation and synchrotron 
x-ray powder diffraction. Substituting Fe by Co weakens the ordered magnetic state through both a 
suppression of Tn and a reduction in the size of the ordered moment. Upon further substitution of 
Fe by Co the high sensitivity of the muon as a local magnetic probe reveals a magnetically disordered 
phase, in which the size of the moment continues to decrease and falls to zero around the same point 
at which the magnetically-driven structural distortion is no longer resolvable. Both magnetism and 
the structural distortion are weakened as the robust superconducting state is established. 

PACS numbers: 74.90.+n, 74.25.Ha, 76.75.-H 



I. INTRODUCTION 

The phase diagrams of iron-based superconductors 
demonstrate that, as for the cuprates, the superconduct- 
ing state generally emerges from a magnetic parent com- 
pound under chemical substitution^— . However, in con- 
trast to the cuprates where the magnetic state is due 
to electrons localized in a Mott insulator, in iron-based 
superconductors it arises from an instability of the delo- 
calised Fe d-band electrons 4 . It is important to under- 
stand how this magnetic state evolves into the supercon- 
ducting one and how the two states compete, coexist or 
mutually exclude one another. 

In this paper we study the "111" arsenide compound 
NaFei-zCo^As using SQUID magnetometry, muon spin 
rotation (/iSR) and synchrotron x-ray powder diffraction 
(XRPD) experiments. Our results allow us to follow in 
detail how the magnetic state collapses as Fe is substi- 
tuted by Co, first by weakening the magnetic long-range 
order (LRO) and structural distortion and then by sup- 
pressing magnetic fluctuations. Our data demonstrate 
the decisive influence of the magnetism on the system and 
show how its disappearance correlates with the strength- 
ening of the superconducting phase. A summary of our 
findings is shown in the phase diagram in Fig. [T] (which 
includes some data points from Ref. 5). It is notable that, 
in comparison with other pnictides, very small amounts 
of cobalt on the Fe site result in a complete suppres- 
sion of the magnetic state. For example, the optimal 
T c is obtained with a substitution of ~ 1.5% Co on 
the Fe site^ in NaFei-jcCo^As, compared to 3 - ~ 6.5% 
Co in Ba(Fei_ a; Co a ;)2As2. This means that significantly 
smaller structural and electronic changes are being made 
to the stoichiometric system thereby reducing any effect 
of inhomogeneous substitution on the properties under 
investigation. 

For this study, a series of NaFei-^Co^As compounds 
were synthesised according to the methods described in 



Ref. [H and details of structural characterisation and pu- 
rity checks can be found there. 
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FIG. 1: Phase diagram for NaFei-ajCo^As compiled using 
data from XRPD (T s ), ZF-^iSR (Tn) and SQUID magnetom- 
etry (T c ) and showing regions of antiferromagnetism (AF), 
superconductivity (SC) and inhomogeneous magnetism (IM), 
as well as their regions of coexistence. The shading of the 
SC region corresponds to the volume fraction estimated from 
susceptibility data. 



II. SQUID MAGNETOMETRY 

Magnetic susceptibility measurements were carried out 
on a Quantum Design MPMS-XL SQUID magnetometer 
under zero-field cooled and field cooled conditions in a 
measuring field of 5 mT. The evolution of zero-field sus- 
ceptibility with temperature for a representative set of 
samples is shown in Fig. [5] (a more complete data set 
can be found in Ref. |5| and the superconducting T c val- 
ues for all studied samples are indicated in Fig. [TJ . Also 
shown is the evolution of the estimated superconducting 



volume fraction with doping [Fig. [2](inset), shown also 
by the shading of the superconducting phase in Fig. [T], 
which suggests that samples with x > 0.0125 are fully 
superconducting. No traces of magnetic impurities, such 
as Fe, were identified by these measurements in any of 
our samples. 
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FIG. 2: Susceptibility data for various NaFei_ x Coa;As sam- 
ples. Inset: estimated volume fractions for all samples with 
< x < 0.02. A full superconducting volume fraction is 
established around x w 0.0125. 



III. ZERO-FIELD pSR 

To probe the magnetic order parameter in 
NaFei-zCozAs, zero-field (ZF) /iSR measurements 
were made on samples with < x < 0.02 using the GPS 
spectrometer at the Swiss Muon Source and the MuSR 
spectrometer at the ISIS facility. Fig. [3] shows example 
spectra measured at T = 1.5 K. For small x, we observe 
oscillations in the time dependence of the muon spin 
polarization (via the positron decay asymmetry), which 
are strongly indicative of magnetic long-range order 
(LRO). These become less distinct with increasing x 
and are not discernible for x > 0.0125. The frequencies 
of the oscillations are proportional to the local magnetic 
field at the muon site(s) and scale with the ordered 
moment, so any one of them may be considered an 
effective magnetic order parameter for the system. 

To extract the order parameter as a function of temper- 
ature and locate Tn, the oscillating asymmetry function 
for the materials with x < 0.0125 were fitted to the ex- 
pression A(t) = Ym=i A* cos(2nUit) exp(— Xtf), with the 
frequencies held in fixed proportion. For all samples, ex- 
cept stoichiometric NaFeAs, it was found that two fre- 
quencies were required to fit the data (i.e. n = 2) with 
one around 20 times larger than the other. In the un- 
doped compound the addition of a third component with 
a frequency very similar in value to the the other low 
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FIG. 3: Comparison of the ZF /iSR asymmetry spectra for 
NaFei-zCozAs, measured around 1.5 K. 



frequency component, slightly improved the fit-. As the 
damping rates, Xi, of all oscillations increase with x, sug- 
gesting a rising level of magnetic disorder, we suggest 
that the third frequency is not resolvable in the doped 
compounds. For all samples with < x < 0.0125, it 
was found that the amplitude of the high frequency com- 
ponent, A\, made up between 55-65% of the total oscil- 
lating asymmetry and the low frequency component (s) 
accounted for the remainder. Muon precession frequen- 
cies for this series, extrapolated to zero-temperature, are 
shown in Table U along with those observed in other iron- 
arsenide systems for comparison. It is notable that the 
high frequency component in NaFeAs is markedly lower 
than that observed in the other iron-arsenides listed, and 
we return to this point later. 



TABLE I: Muon precession frequencies extrapolated to T — 
K, along with amplitudes, for various FeAs-based com- 
pounds (* indicates majority component). 



Compound 



Frequencies (MHz) 



Ref. 



NaFeAs 

LaFeAsO 
SrFeAsF 

BaFe 2 As 2 
SrFe 2 As 2 

FeAs 



10.9*, 0.9, 0.5 

23*, 3 
22.2* , 2.0 

28.8*, 7 
44* , 13 

38.2, 22.7* 
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[6] 
17] 

18] 
[9] 

[10] 



For the remaining samples (0.015 < x < 0.02) a fast- 
relaxing component is observed at early times that disap- 
pears as x increases. These data are best fitted to a relax- 
ation function A(t) = A s i ow exp(—a 2 t 2 )+Af ast exp(— Aft), 
which combines a slowly-relaxing Gaussian function, at- 
tributable to the magnetism of static and disordered 
moments probably of nuclear origin (a Gaussian Kubo- 
Toyabe function with the same magnetic field distribu- 
tion fits just as well in this regime), added to a fast ex- 
ponential relaxation. This fast relaxing component is 
most likely attributable to magnetic fluctuations of elec- 
tron spins which quickly dephase muons. The fact that 
these two contributions are summed suggests that only a 
fraction of the muons are dephased by the electronic fluc- 
tuations. This implies that, for x > 0.0125, samples show 
some degree of inhomogeneous magnetism (IM). For the 
materials in this region we observed a decrease in the rel- 
ative amplitude of the fast component, v4f ast , indicating 
that the extent of the magnetically fluctuating regions 
decreases as x increases. Our muon data show no sign of 
signals from the most common impurities in iron arsenide 
compounds^. 
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FIG. 4: (a) Precession frequency (proportional to ordered 
moment) as a function of Neel temperature shown by filled 
circles. In the IM state, the effective ^i(0) is estimated from 
the magnetic contribution to the rms field width and Tn is 
taken to be the temperature at which this contribution sets 
in [see Fig. [7]. Inset: The probability distribution for the 
ordered moment for stoichiometric NaFeAs given the observed 
precession frequency, (b) The estimated ordered moment as 
a function of x in the AF state (filled circles) and IM state 
(open circles). 

The low values of the ZF muon precession frequen- 
cies in NaFeAs, compared to those seen in other iron- 
based superconductors (see Table JJ, indicate that the 
moment is correspondingly low. Neutron scattering mea- 
surements indeed indicate a low moment, 0.09(4) /ib, on 
the Fe site in pureii NaFeAs. Our ZF-//SR data can be 
used to make an independent estimate of this moment, 
along with its variation with x. We can obtain the prob- 
ability P(n-F C \{vi}) [shown in the inset of Fig.@Ja)], that 
is the probability of /^Fc taking a particular value given 
that we observe a set of precession frequencies, {vi}, us- 
ing Bayes' theorem to invert the probability of observ- 
ing a set of precession frequencies given a particular Fe 



moment (which is easily calculable from dipole-field sim- 
ulations). This method (see Refs. [l2l and [ill ) allows us 
to provide a quantitative estimate of /iF C = 0.10(3)/xb, 
assuming the magnetic structure to be as determined in 
Ref.Hl 

ignoring any hyperfine coupling and constraining 
the muon location to be at least 1 A from an atom site. 
This estimate is in good agreement with the value from 
neutron scattering^. The dependence of ^i(0) as a func- 
tion of Tn is shown by the solid circles in Fig. Eta) and 
hence the dependence of /iFc as a function of x is shown 
in Fig. QJb). It is noticeable that ^i(0) (and hence jUFe) 
tends to zero more quickly than Tn. This suggests it is 
the collapse of the moment that destroys the magnetic 
state. 



IV. X-RAY POWDER DIFFRACTION (XRPD) 

To probe the structural distortion in NaFei-^Co^As 
we performed XRPD studies of samples with < x < 
0.025 on the ID31 beamline at ESRF. In Ref. i the dis- 
tortion was characterised as a transition from a high tem- 
perature tetragonal (PA/nmm) phase to a low tempera- 
ture orthorhombic (Cmma) phase, with a ort h = et- 
In order to analyse the variation in the size of this dis- 
tortion with temperature and doping, we fitted data at 
all temperatures to a Cmma phase. In the TOPAS Aca- 
demic software^, the lattice parameters were defined as 
a = \/2atot+^/2 and b = v^atct — S/2, with both a tc t and 
6 allowed to vary. The data were then fitted to this model 
using a structure-independent Pawley refinement^ and 
a typical refinement is shown in Fig. [5] At high temper- 
atures and/or dopings, where S was sufficiently low as to 
only cause peaks to broaden rather than split, a set of 
refinements with S fixed at zero (i.e. using a tetragonal 
model) was carried out. Comparing the statistics of the 
two models allowed us to confirm the temperatures at 
which the structural distortion was no longer resolvable. 
Plots of 5 vs T for all samples with a resolvable distortion 
are shown in Fig. HJa). These data sets have been fitted 
to the function 5 = [1 - (T/T S ) Q ]' 3 to locate T s , with a 
allowed to vary and j3 fixed around the mean-field value 
of 1/2. 

Combining Tn from /iSR, T s from XRPD and mea- 
surements of superconducting T c and volume fraction 
from SQUID magnetometry, allows us to produce the 
phase diagram in Fig. [TJ showing that the superconduc- 
tivity, magnetism and structural distortion coexist for 
x < 0.0125, with the superconducting volume fraction 
increasing with increasing x. The results of our ZF-/iSR 
experiments indicate that the magnetism weakens as x 
is increased, with the system unable to sustain magnetic 
LRO for x > 0.0125. However, IM and the structural dis- 
tortion coexist with bulk superconductivity in the region 
0.015 < x < 0.020. 

Further insight into the effect of the reduction in the 
strength of the magnetism may be obtained by comparing 
the magnetic and structural order parameters, fi(T) and 
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FIG. 5: A Pawley refinement for NaFeo.9sCoo.02 As at 5K. 
Inset: A comparison of the fits to the 022 and 202 peaks in 
Cmma at 5K and 10 K. These peaks converge into the 112 
peak in PA/nmm. Note that A, the splitting in the 29 values 
of the two peaks, is smaller at 5 K than at 10 K, suggesting a 
suppression of the structural distortion at low temperatures. 
This is in agreement with the behavior of 5(T) as obtained 
from fitting the full pattern. 
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S(T) respectively, as shown in Fig. |6] A reduction in 
both ^i(0) and <5(0) is apparent as x is increased, along 
with a marked suppression of these order parameters with 
decreasing T in the more highly doped samples (indicated 
by arrows in Fig. [5]). The suppression of the structural 
distortion is evident in Fig. \5linset). In this figure the 
difference in the fitted 29 values (A) of the orthorhombic 
022 and 202 reflections (into which the tetragonal 122 
peak splits) falls on cooling from 10 K to 5 K (having 
increased on cooling from T s to 10 K). We note that the 
suppression of both S and v at low T is observed in those 
samples with a full superconducting volume fraction, in 
agreement with what is seen in other pnictide families*™*. 



V. TRANSVERSE-FIELD ^SR 

To further probe the region above x — 0.015 we 
have used transverse field /xSR, which also provides a 
method of probing the superconducting state of type- 
II superconductors^. A magnetic field is applied per- 
pendicular to the initial muon-spin direction, which also 
produces a vortex lattice in the material's superconduct- 
ing phase. This results in a distribution in the local 
field across the sample, whose rms width B rnls can be 
extracted from the muon precession. Contributions to 
B rms arise from superconductivity, static electronic mag- 
netism and nuclear moments (the latter being small and 
temperature-independent) and add in quadrature. The 
measured temperature dependence of B lms is plotted in 
Fig. [7] and for each sample studied we find that £? rms 
rises sharply on cooling through the superconducting T c 
(found independently from SQUID magnetometry) con- 



FIG. 6: The evolution of (a) structural and (b) magnetic 
order parameters with T and x. S is defined as the difference 
between the a and b lattice parameters in the low-temperature 
Cmma phase and v\ is the largest observed muon precession 
frequency in zero- field. The blue line sections show where 
the sample superconducts and the shade of blue indicates the 
volume fraction obtained from SQUID magnetometry. The 
arrows indicate the suppression of S and v\ in the more highly 
substituted samples. 

sistent with a contribution to the broadening from the 
vortex lattice. In addition, we observe a contribution 
to -B rms from magnetism which decreases steadily as x 
increases from 0.015. For x — 0.025 the magnetic con- 
tribution is negligible and the temperature dependence 
of -B rms is consistent with that expected for conventional 
superconducting order [solid line in Fig. [TJd)] . It is likely 
that in the materials for 0.015 < x < 0.025 the IM is as- 
sociated with a non-zero moment on the Fe site, but the 
moment is insufficient to stabilize a fully ordered mag- 
netic state. We then associate the increase of -B r ms just 
above the superconducting T c with the onset of the IM 
state. We also observe a peak in B lms at around 17 K 
(most noticeable in the sample with x — 0.015 but weak- 
ened and broadened for higher x) which may be due to 
increased correlation between the static moments (but 
which stops short of LRO, as we do not observe a pre- 
cession signal in ZF-/1SR). This peak is also observed in 
the ZF amplitude of the fast-relaxing component for the 
x = 0.015 sample [Fig. 0a)], lending further weight to 
this interpretation. 

We can extract an estimate of the magnetic contri- 
bution to B ms in the IM regime by assuming that a 



magnetic contribution and a non-magnetic contribution 
(the latter assumed to be that of the x — 0.025 sample) 
add in quadrature. This rms field width can then be 
used to provide an estimate of the precession frequency 
[(7 M /27r)i? rms ] that would have been observed if the mo- 
ments in the IM state ordered, and we plot these as open 
circles in Fig. QJa). By assuming the same scaling be- 
tween frequency and ^Fc as for the AF-ordered region, 
we can use these values to estimate /zf c which is plotted 
in Fig. IDJb) using open circles. These data show that 
as x increases, the moment on the Fe site undergoes a 
process of steady suppression, initially resulting in a loss 
of AF order as the order loses long-range coherence but 
then collapsing further in the IM state. 
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FIG. 7: The width B rms vs T for samples close to the mag- 
netic LRO region. For x — 0.0150 (a) the unusual peak in 
B rms matches one observed in the zero-field fast relaxation 
amplitude, A{ ast . The grey line indicates T c from SQUID 
magnetometry (see Fig. [TJ. The black line in (d) is a fit to a 
conventional SC order parameter temperature dependence. 



VI. DISCUSSION 

It has recently been suggested that the structural 
transition in the iron pnictides is driven by magnetic 
fluctuations^ 8 -, rather than by magnetic order. The 
fluctuations could be of different types: spin-nematic 
fluctuation o 18 ' 19 , ferro-orbital fluctuations or critical 
fluctuation o 22 ' 23 . Our discussion will focus on the spin- 
nematic model of Ref. [IH which we believe explains the 
main features of our data. The theoretical picture is mo- 



tivated by the magnetic structure found in these materi- 
als, where Fe spins are coupled antiferromagnetically for 
Fe-Fe bonds along one diagonal of the two-Fe tetragonal 
unit cell and ferromagnetically along the other diagonal. 
This results in two coupled antiferromagnetic sublattices 
with coplanar staggered magnetization. The sublattice 
coupling is related to a Z2 symmetry, which allows two 
possible orientations of stripes in the magnetic structure. 
This symmetry is broken by nematic order, which can be 
present with zero sublattice magnetization and whose on- 
set may therefore occur at temperatures above any mag- 
netic ordering temperature. The result of this nematic 
order is to bias the magnetic fluctuations along a par- 
ticular stripe orientation. Crucially, the nematic order 
parameter couples to the lattice, causing bonds between 
neighboring (anti)parallel spins to (expand) contract and 
this leads to the observed structural distortion. The pre- 
requisite for nematic order is a large, but finite, mag- 
netic correlation length, whose occurrence at elevated 
temperatures is probable for materials in both the mag- 
netically ordered and IM region of the phase diagram of 
NaFei-zCozAs. 

From this point of view it is magnetic fluctuations, bi- 
ased by nematic order, that not only lead to magnetic 
order, but also drive the structural distortion. Substi- 
tuting Fe by Co causes a reduction in the strength of 
these magnetic fluctuations and this suppresses both the 
magnetic order and the structural distortion. Our mea- 
surements show that on increasing x the structural and 
magnetic order parameters weaken until the fluctuations 
can no longer sustain magnetic order above x > 0.0125. 
Further substitution of Fe by Co continues to weaken the 
fluctuations and this coincides with the weakening of the 
structural distortion order parameter until fluctuations 
can no longer sustain the structural order parameter for 
x > 0.0225. This weakening of the magnetic fluctuations 
is coincident with a strengthening of superconductivity in 
NaFei-zCo^As. This is also consistent with the nematic 
order picture, which predicts a competition between su- 
perconductivity and magnetism. Specifically, the onset 
of superconductivity is predicted to lead to a reduction 
in the static part of the magnetic susceptibility 18 which 
weakens the magnetic order and spin correlations. This is 
seen in our data, where we observe the low T suppression 
of (i) the magnetic order parameter in the x = 0.0125 
material (where long-range magnetic order is only just 
sustained) and (ii) the structural order parameter for 
x > 0.15 (where the structural order parameter has been 
weakened). We note that these samples have close to 
full superconducting volume fraction and the dip in the 
order parameter is only observed for T <C T c , where the 
superconducting order parameter has become sufficiently 
strong. 

It is interesting to compare these results with those ob- 
tained on another "111" superconductor, LiFeAs 2 ^— . A 
result of accommodating the smaller Li + ion (rather than 
Na + ) between the FeAs layers is that the edge-staring 
FeAs4 tetrahedra in LiFeAs are very compressed in the 



basal plane. Superconductivity occurs in the stoichio- 
metric material LiFeAs, in contrast to other iron arsenide 
superconductors for which doping away from the formal 
oxidation state of Fe 2+ or the application of hydrostatic 
pressure is required to induce superconductivity. Fur- 
thermore, stoichiometric LiFeAs does not appear to show 
static magnetism^! (in contrast to NaFeAs) . The present 
study demonstrates that substituing a small amount of 
Co on the Fe site strengthens the superconducting state 
in NaFeAs, but in LiFeAs substitution of Fe by small 
amounts of Co or Ni results in a steady lowering of the 
superconducting transition temperature^. For LiFeAs, 
T c is lowered monotonically at a rate of 10 K per 0.1 elec- 
trons added per formula unit irrespective of whether the 
dopant is Co and Ni, and at higher substitution levels 
superconductivity is completely suppressed. (Co and Ni 
have a similar structural effect as a function of the level of 
doping, but Ni adds twice as many electrons as Co.) The 
number of electrons added per formula unit is also the 
determinative quantity for NaFeAs, but here the super- 
conducting state is first strengthened (as the magnetic 
state is destroyed) and then weakened (see Fig. 1). 

Finally, the behavior of the supcrfluid stiffness as a 
function of T c in LiFeAs derivatives is markedly different 
from that of other pnictides^, including the isostructural 
NaFeAs-derivatives considered here and the "1111" and 
"122" classes^. Most of the pnictides exhibit behav- 
ior more similar to the hole-doped cuprates, but LiFeAs 
derivatives resemble much more closely the electron- 
doped cuprates. The differing behavior of LiFeAs may 
result from the underlying structural difference intro- 



duced by the small Li + ion, resulting in a band structure 
which does not favor Fermi surface nesting^ so that the 
magnetic instability found in the other isoelectronic (i.e. 
undoped) iron arsenides does not compete successfully 
against superconductivity in LiFeAs. 



VII. CONCLUSION 

In conclusion, we have presented a picture of 
NaFei-zCo^As in which a magnetic interaction drives 
both magnetic LRO and a structural distortion and 
weakens with the subsitution of Fe by Co, causing a de- 
crease in the size of the ordered magnetic moment. Mag- 
netic LRO becomes unsustainable above x — 0.0125 but 
the structural distortion persists until the moment in the 
IM state is no longer strong enough sustain it. The IM 
state may be characterized by a long, but finite magnetic 
correlation length and nematic order, which biases the 
magnetic fluctuations along stripes. Competing directly 
with magnetism is superconductivity which, on reaching 
full volume fraction, weakens the magnetism dramati- 
cally and is seen directly to depress the magnetic and 
structural order parameters at low temperature. 
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